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We relate the MSW effect to the efective absorption of the 
electronic collective motion energy by retaining the imaginary 
part of the index of refraction associated with the charged- 
current scattering and show that the small angle MSW so- 
lution to the solar neutrino anomaly can be used as a probe 
of the physical conditions of the solar interior if it is correct. 
We find that the constraint on the absorption imposed by 
the small angle MSW solution and the theoretical estimate 
of the absorption by the Boltzmann kinetic theory are consis- 
tent, which shows that a consistent theoretical picture can be 
developed when plasma absorption processes are taken into 
account. 
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If, based on its apparent success, the small angle MSW 
solution is accept as the solution for the solar neu- 
trino problem , it can be used as a probe of physical 
conditions in the solar interior. Since the MSW effect ]Tc| ] 
can be attributed to the resonant scattering of neutrinos 
off electrons in the solar interior, it must be sensitive to 
the absorption of not only the neutrino energy but also 
the electron energy. The effective absorption can be rep- 
resented in terms of complex indices of refraction. In 
order to take into account the absorption, we rewrite the 
index of refraction n e associated with the charged-current 
scattering in such a suggestive form pi : 



G, 



(1) 



where Gp, a — e 2 , A c = m" 1 , E and u e are, re- 
spectively, the Fermi constant, fine structure constant, 
electron Compton wavelength, neutrino energy and elec- 
tronic plasmon energy, noticing that we have set K = c = 
1. Since the MSW effect depends on an effective density- 
dependent contribution to the neutrino mass, longitudi- 
nal plasmons or Langmuir waves due to electrostatic os- 
cillations of free electrons in the solar plasma make sense 
lM3: 



= « + 3y T 2 e fc' 2 ) 1/2 + *r', 



(2) 



where u pe = (4irN e e 2 /m e ) 1 ' 2 , V Te = {T/m e ) 1 ' 2 , k' 
and r' are, respectively, the electronic plasma frequency, 
characteristic thermal velocity of electrons, plasmon mo- 
mentum (or wavenumber) and effective absorption co- 
efficient of the plasmon. The fact that uj e is complex 
means that the plasmon or the coherent motion of elec- 
trons in the plasma e _l " c * will decay. If the plasma is 



not uniform, all dependent variables will depend on the 
location. Obviously, if k' — and V = 0, then uj e — lu pb , 
and hence E(n e — I) = G F N e , as assumed in the MSW 
theory |io|Jl]-|^|. The purpose of this letter is to investi- 
gate how the small angle MSW solution constrains V in 
the solar plasma. 

For simplicity, we consider two neutrino flavors. If 
we use the flavor eigenstates \v e ) and \v p ) as the basis, 
the time evolution of the neutrino state vector \v(x)) — 
a e {x)\v e ) + a M (x)|f M ) in matter in the relativistic limit is 
governed by the equation (h]|| 



dx 



Ay / M{x) sin 26V 
2 I sin 26V ~M(x) 



(3) 



where x = ct, M(x) = y/2E(n e (x) - 1)/Ay - cos 26V, 
Sm 2 = m|— m 2 , Ay = Sm 2 /2E, By is the vacuum mixing 
angle. M r {x) = ReM(x) and Mi(x) = ImM(x) can be 
approximately cast as follows: 

M r (x) w V2G F N e (x){l + 3fc 2 - r(x) 2 )/Ay 

-cos 26V, (4a) 
Mi(x) w 2V2G F 7V e (a;)r(2;)/Ay, (4b) 

where k — k'/kn, r(x) = V [x) /to pe {x), and kn = 
Upe/Vre- We have assumed fc 2 <C 1, which is reason- 
able because when fc ^ 0.3 the Landau damping |ll] |TJ] 
will cut in. The fact that M(x) is complex implies that 
there is no N e (x c ) so that M(x c ) = no matter N e (x) is 
large or small. Nevertheless, if the absorption T is small 
enough, the level crossing still occurs when 



M r (x c ) > 0. 



(5) 



Obviously, k favors resonance, while T disfavors reso- 
nance as it should. 

At a free electron number density, N e (x), the light (L) 
and heavy (H) local mass eigenstates [[15|,[§ are 



\vl{x)) = cos9(x)\v e ) ~ bwl9(x)\v ij ,) 
\vh{x)) — sin0{x)\v e ) + cos0(x)|f /x ), 

which have complex eigenvalues ±iA(x), where 

A(x) = A v [M(x) 2 + sin 2 26V] 1 / 2 , 
and 9(x) satisfies 

sin 26V 



sm26(x) = 
cos26'(a;) = 



A(x)/Ay' 
-M{x) 

A(x)/Ay' 



(6a) 
(6b) 



(7) 

(8a) 
(8b) 
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Since 



A(x c ) = A y (sin 2 29y - 4r 2 cos 2 29 v ) 1/2 , 



(9) 



the absorption will affect the width of the resonance: 
2A(x c ), where T c — T(x c ). The corresponding resonance 
distance [|| is 



5x = 2 



Ne 



(tan 2 29 v - 4E 



2x1/2 



(10) 



where N e (x) = dN e (x) / dx is the electronic density gradi- 
ent. Because the small angle MSW solution gives the best 
fit to existing solar neutrino data, sin 2 29y ~ 8 x 10~ 3 , 
and Sx > 0, we obtain 



ir,. 



< itan26>y 



0.04, 



(11) 



which is the upper limit of |r c | in the solar interior where 
the resonance takes place. 

The fact that the eigenvalues are complex implies that 
the neutrino states will decay, too. In order to see 
this we express the neutrino state vector \v{x)) in the 
matter in terms of the local mass eigenstates \v(x)) — 
aii{x) \vh(x)) + (Il{%)\v Consequently, the evolu- 
tion equation becomes 



dx 
where 



|A(a;) ia(x) 
—ia{x) — ^A(x) 



(12) 



a(x) 



A v V2G F N e (x) sin 29\ 
~2 A(a;) 2 



enforces mixing of the mass eigenstates governed by the 
density gradient, where we have assumed dT /dx = 0. 
When the off-diagonal elements can be neglected with 
respect to the diagonal elements, i.e., when 



j(x) 



A(x) 



a(x) 
sin 2 29 



\A(x) 



cos26>y A 2 / sin 3 26»v 



N e (x) 



N r . 



» 1, 



the neutrino will propagate adiabatically through the 
matter. The absorption tends to weaken this condition, 
for example, near the crossing point, it reads 



7< 



= 7c °(l-4r 2 cot 2 2M 3/2 <7 c °- 



(13) 



where 7° = 7(x c ;r c = 0). Since the detected solar 7 Be 
neutrino flux |5|-[)| is much less than the predicted B , the 
7 Be neutrinos must propagate adiabatically in the sun, 
the absorption correction factor thus should not depart 
from unity significantly, which demands 



ITJ « 0.04. 



(14) 



Of cause, the absorption also affects the adiabatic bound- 
ary pl| of 7 C ~ 1 in the 5m 2 / E — sin 2 29y plane. 
When 7 C ^> 1, the adiabatic states are 



\H) 



e- ia{xo '' x)+b{xo '' x) \v L {x)), 
e ia{xo ' x} - b(xo '' x) \is H (x)), 



where 



a(xo; x) = i / A r {x')dx 



b(x ; x) = i / Ai(x')dx\ 



(15a) 
(15b) 



(16a) 
(16b) 



where Xq is the production point of the basis states, 
A r (a;) = ReA(x), Ai(x) = lmA(x). Obviously, they 
will decay provided that b(x) ^ 0. If b(x) > 0, the local 
heavy mass eigenstate will continuously jump down into 
the local light mass eigenstate, while the local light mass 
eigenstate will continuously jump up into the local heavy 
mass eigenstate; if b(x) < the processes are reversed. 
Expressions (0), (pq), ( |16b| ) and condition ( pT[ ) show 
that b(x) < if both the plasmons are indeed absorbed 
(T < 0) and the plasma density is supercritical (M r > 0) 
or both the plasmons are excited (T > 0) and the plasma 
density is subcritical (M r < 0). However, it is still conve- 
nient to use these states as the basis states in the region 
for which there are no transitions Under the adi- 
abatic approximation, we have \v{x)) = a\\L) + a2[H) 
in which the linear combination coefficients ai (i=l, 2) 
are determined by the initial condition, hence the aver- 
age probability of detecting an electron neutrino at the 
Earth with the initial condition \v{xq)) — \v e ) is 



K e d = |(1 + cos2(9 cos2(9y-) cosh26 v 
+ i(cos 26»o + cos2^v)sinh 2b v , 



(17) 



where 9q = 9(xq; T = 0), by = b{x^\ R®). This expres- 
sion will reduce to the well-known adiabatic probability 
if r = or by = 0, noticing that sinh26y = and 
cosh26y = 1 when by — 0. This shows that the absorp- 
tion also affects the detection probability. 

However, it is necessary to go beyond the adiabatic 
approximation because the solar pp and 8 B neutrinos are 
suppressed partially We follow Parke j| to do so. 

Eqs. (|l0|), ([l3]) and ( |l7| ) show that the critical region for 
nonadiabatic behavior occurs in a narrow region [x_ , x+] 
(for small 9y) surrounding the crossing point x c and that 
this behavior is controlled by both the density gradient 
[N e /N e ] Xi . and the absorption T c at the crossing point, 
where x± = x c ±j5x. See Fig. 1 of Bethe || for graphical 
understanding. When an initial electron neutrino at x$ 
approaches to the inner boundary of the nonadiabatic 
region, its state becomes 



\v{x^)) = cos# e 1 
+ sin 9 e 



l - +h -\v L {x-)) 
a -»-\v H (x-)), 



(18) 
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where a_ = a(xo] X—), 6_ = b{x^\ xJ). As the neu- 
trino goes through the nonadiabatic region to reach its 
outer boundary x+, we have the mixed states according 
to Parke § 

\v L {xJ)) -> a 1 \v L (x + )) + a 2 \yH{x+)), (19a) 
\v H {x-)) -» -a* 2 \u L (x+)) + a\\v H {x + )), (19b) 

where a, (i=l, 2) are determined by the nature of the 
transition point and satisfy |ai| 2 + | a-2 1 2 = 1. From 
the upper boundary on, the eigenstates \vl{x)) and 
\vh(x)) will evolve adiabatically starting with \vl(x + )) 
and \vh{x+)) respectively. Therefore, the state vector of 
the neutrino reads as follows: 

\v{x))=A\v L {x))+B\v H {x)), (20) 

in the detection region x > x c , where 

A = ai cos0 o e~* A++B+ - a* 2 sin6 e~' A - +B - (21a) 
B = a 2 cos9 e lA -- B - + a[ sm9 e lA +- B + (21b) 

in which A± = a+ ± a_, B± = b + ± 6_, a + = a(a;+; x) 
and 6+ = 6(a;+; x). 

Substituting Eqs. ( |6a| ) and ( |6b| ) into Eq. (|20|), one can 
find the amplitude for producing, in the solar core xq , and 
detecting, on the Earth, an electron neutrino after pas- 
sage through resonance: A e = AcosOy + BsinOy- Thus 
the probability of detecting this neutrino as an electron 
neutrino after averaging over both the production and 
the detection positions is given by 

P Ve = |ai| 2 (cos 2 O cos 2 9 v e 2B+ + sin 2 6 sin 2 6Ve" 2B +) 
+ |a 2 | 2 (cos 2 6>osm 2 6Ve 2S - + sin 2 9 a cos 2 9 v e- 2B - ), 

(22) 

where b + = b(x+; Rq)- If B± = 0, we reproduce the 
well-known Parke formula : 

pParke = 1 + ( 1 _ j cos 2 q q ^ ^ ^3) 

where |ai| 2 = 1 — 1 0,2 1 2 and P x = |a,2 1 2 is the probability of 
transition from vh(x-)) to (or vice versa) as the 

neutrino goes through the nonadiabatic region. Parke H 
and Haxton |lq| have worked out P x by using a linear 
density profile, it is natural to generalize it to our case 
with T ^ 0: 

P*=exp(-| 7c ). (24) 

Because the small angle MSW solution with B± = 
gives the best fit to existing solar neutrino data, the solar 
neutrino experiments require \B±\ -C 1. Using the mean 
density pq ~ 1 g cm -3 , the mean number of electrons per 
nucleon Y e = |, we can estimate \B±\ ~ 10 2 |r|, where 
r is the mean absorption coefficient in the solar interior. 
Consequently, the small angle MSW solution demands 



|r| < icr 2 . (25) 

If we assume that the resonance region is near the produc- 
tion region, then b- w 0, hence B + ps B_ w b + . Since 
Aj(x) fa ~r(x)2V2G F N e (x)cos26'(x), where 6'{x) = 
9(x; r = 0), and cos 29' (x) > when x > x+, then 
Ai(x) > when T(x) < 0. Therefore, B± > when 
r(x) < 0, but _B± < when T(x) > 0. Since numerical 
evaluation shows that Eqs. (^) and (|23| ) give the same re- 
sults when B± < 10 -3 , we may infer that T(x > x c ) < 
and |r| < 10 -5 in the solar interior. 

The Boltzmann kinetic theory formulated in the Boltz- 
mann equation jllj can be used to estimate the effective 
absorption coefficient of the electronic energy. The effec- 
tive energy transition rate Z e i and momentum transition 
rate R e i from electrons to ions per electron in the elastic 
scattering process p e + p; ^ p' e + p- (Chapter 3 of Jil"|) 
can be estimated through Tayler-expanding these rates 
near local thermodynamic equilibrium between electrons 
and ions as follows 

Zei = -i4(T e -Ti), (26a) 
R« = -v%(T e - T t )(m e u e - m lUl ), (26b) 

where and V = — z/g(T e — Ti) are the effective en- 
ergy and momentum transition (or absorption) coeffi- 
cients of electrons in the plasma, u e /j) is the coherent 
motion velocity of electrons (ions) , noticing that Z^ = 
and = when electrons and ions have the same 
temperature T e = Ti. Since the effective energy and mo- 
mentum transition rates are dependent on each other, 
r = -v e z l (T e -Ti)/T e . Therefore, we need only the effec- 
tive energy transition coefficient , which may be esti- 
mated by the thermal timescale of the sun, the timescale 
that the stored thermal energy of the sun is used up via 
radiation at the present radiation power. The thermal 
timescale is equal to the gravitational timescale, which 
is the ratio of the gravitational energy to the total lumi- 
nosity f§ 

tgra^ty ~ GM^/RqIq « 10 7 yr. (27) 

So we may have 

r'~ -lO-^Te-TiVTosec- 1 . (28) 

Since u> pe 3> 1 sec -1 and T e > Ti in the solar interior, 
we know T = T'/oj pe is negtive and |T| <C 10~ 5 in the 
solar interior. This shows that the small angle MSW so- 
lution to the solar neutrino problem and the Boltzmann 
kinetic theory are consistent. A treatment of the plasma 
processes including absorption thus verifies that a con- 
sistent theoretical picture can be developed. 
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